Normal Modes of the B=4 Skyrme Soliton 
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The Skyrme model of nuclear physics requires quantisation if it is to match observed nuclear 
properties. A simple technique is used to find the normal mode spectrum of the baryon number 5 = 4 
Skyrme soliton, representing the a particle. We find sixteen vibrational modes and classify them 
under the cubic symmetry group Oh of the static solution. The spectrum possesses a remarkable 
structure, with the lowest energy modes lying in those representations expected from an approximate 
correspondence between Skyrmions and BPS monopoles. The next mode up is the 'breather', and 
above that are higher multipole breathing modes. 



I. INTRODUCTION 

In the Skyrme model of nuclear physics, both pions 
and nucleons are represented by a single scalar SU(2) 
group valued field, U(x). Pions occur as field quanta, 
while baryons are instead represented as topological soli- 
tons. The classical Skyrme theory, with a simple quanti- 
sation of the spin and isospin collective modes, provides 
a description of nucleons and the A resonance in modest 
agreement with experiment 

Applying the Skyrme model to larger nuclei and to 
nuclear matter is an even more interesting proposition, 
since with no additional free parameters one could com- 
pare the theory with the binding energies and gamma ray 
spectra of all nuclei. There has been progress in under- 
standing the structure of Skyrme multisolitons [Q, 
but it is clear that unless quantum fiuctuations about 
the static solutions are included, there is little chance of 
success. For real nuclei the relative kinetic energies of 
the nucleons in the ground state are large, so a quantisa- 
tion at least of a number of degrees of freedom equal to 
the number of nucleon coordinates is essential. Note that 
the simple collective coordinate quantisations of spin and 
isospin include effects of order , while ignoring ef- 
fects of order fi. 

Recently there has been some progress in quantisation, 
based upon Manton's notion of representing low energy 
solitonic excitations as motion on a finite dimensional 
space of moduli. A study of the deuteron by Leese et. al. 
1^ using an instanton approximation for the field config- 
urations, gave encouraging agreement with experimental 
properties, but only included two of the four expected 
vibrational modes of the deuteron. Walet |^ extended 
this treatment to estimate all vibrational frequencies for 
the deuteron and triton, again using the instanton ap- 
proximation. 

In this Letter we follow a different track. We di- 
rectly compute the low energy normal modes, finding 
their frequencies and the representations they lie in un- 
der the static soliton's symmetry group. The frequencies 
and representations provide a coordinate-independent 



description of the configuration space around the static 
solution; in the harmonic approximation they determine 
the quantum vibrational energy levels. Our results pro- 
vide new insight to the moduli space approach, since the 
representations for the lowest frequency modes turn out 
to be just those expected from a recently understood ap- 
proximate correspondence between Skyrmions and BPS 
monopoles Furthermore, the success of these calcula- 
tions allows one to contemplate going beyond the moduli 
space approximation and performing a full semiclassical 
quantisation of the field theory. This is an attractive goal, 
since the Skyrme theory could then be incorporated in 
the framework of chiral effective Lagrangians jl^ , allow- 
ing a unified treatment of mesons, baryons and higher 
nuclei. 



II. METHOD 

Since the SU(2) manifold is a 3-sphere, we represent it 
in terms of a scalar field (/) € R^, with (Z)*^;/)" = 1. In 
terms of this field, the Skyrme Lagrangian density is 

C = \d^(t) ■ d^'cl) + ojl(t>^ + \{(t> ■ 4) - \) ^ 



with A a Lagrange multiplier field. Here, length and 
time are in units of -4^, energy in units of In 
these rescaled units, the only remaining parameter is 
ti)^ = the oscillation frequency for the homogenous 
pion field. For the most part we have adopted the 'stan- 
dard' value Q of ojtt .526, although we have also per- 
formed calculations at twice this value. 

The mixed space and time derivative terms in the 
Skyrme Lagrangian make numerical solution in general 
difficult. Isolating the time derivative terms one has 

C^^r K'^'id.cl,) / ~V{cf>, d.cl,) + X{cj, . cf, - 1) , (2) 

where K"'' = S^^l + ((9^0)^) - d^cj)" dicj)'' is a local inertia 
matrix, and V{(j), dicj)) is the potential. In general, iC"* 
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is time dependent, but for small perturbations around a 
static solution 4>st{x) we write 0(a;, t) = 0st(a^) + e{x, t), 
e 1 and to second order in e the Lagrangian is 

£ = 1 ^afc (5^^^) _y^^^ Q^^^ + A(0 . ^ - 1) (3) 

This Lagrangian leads to classical field equations: 



dV 



(4) 



where the matrix K°-'^{x) is taken to be its value at the 
static classical solution. Equation (4) closely approxi- 
mates the Skyrme equations for fields near a static clas- 
sical solution, precisely the desired regime for studying 
soliton normal modes. 

In order to numerically solve the field equations, we 
discretise the action, (||), using a diagonal differencing 
scheme for the four spatial derivative terms, achieving a 
high degree of locality and second order accuracy in both 
spatial and time steps. The numerical code conserves en- 
ergy and baryon number to within 1 part in 10^ over the 
course of extremely long (50k timestep) runs. Periodic 
boundary conditions are used. We first create the appro- 
priate minimal energy static solution by straightforward 
time evolution from four-Skyrmion initial conditions. A 
simple relaxation procedure sets the field momenta zero 
each time the kinetic energy reaches a maximum. We 
find the fields rapidly converge on the minimum energy 
configuration. K°'^ is set equal to 5°"^ in this part of the 
calculation, since this does not affect the final static so- 
lution. 

Next we slightly perturb the fields and evolve them 
forward again, but now using the full inertia matrix 
K^'^idcji^t)- The evolving field is 

(j){x,t) = (j)st{x) + ^ e„(5„(a?)cos(w„i) + 0(6^) (5) 

modes 

where the functions Sn{x) G R^, obeying 5n{x) ■ (j)st{x) = 
0, are the normal modes, each excited with amplitude e„. 
The normal mode frequencies a;„ are found by Fourier 
transforming (l){x, t) with respect to time at any point in 
the box, and plotting the resulting power spectrum. 
The space of perturbations has a useful inner product 
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S1{x)K''''{x)S^{x) d^ 



(6) 



box 



which is zero for normal modes 6i and S2 if wi lu2. The 
inner product allows one to determine the degeneracies of 
the normal mode frequencies and the representations in 
which the modes transform under the soliton's symmetry 
group. 
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FIG. 1. Fourier power spectrum for perturbations around 
a B=4 soliton. Frequency is in Skyrme units, power scale 
is arbitrary. Note that the frequency of homogeneous pion 
oscillations is o;^ = .526 in these units. 



III. RESULTS FOR THE B=4 SOLITON 

We have applied this technique to the case of the 
Skyrme soliton with baryon number four: the a parti- 
cle. The a particle provides an especially simple case 
for quantisation, because the ground state possesses zero 
angular momentum and isospin. The static soliton has 
cubic symmetry; its energy and baryon number density 
concentrate along the edges of a cube The full 48 
dimensional cubic group of symmetries Oh (for notation 
see |ll|) is generated by 90 degree and 120 degree rota- 
tions, and parity /. After an appropriate global isospin 
rotation, the action of these group elements on spatial 
coordinates and pion fields ((/)° = ("'7^)) is as follows [0: 
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{TT^,7T'^,7T^){x,y,z) (tT^ , TT^ , TT^ ) (y, Z, 
{TT^,7T'^,7T^){x,y,z) (tt'^ , Tt"^ , TT^){~X , ~y, - z) 



with TT^ — \{t^^ — 27r^ — 27r^) etc. From this it is straight- 
forward to check that a homogeneous pion field falls into 
the two dimensional representation and the one di- 
mensional representation A^, where superscripts indi- 
cate parity. 

A typical power spectrum of the perturbations is shown 
in Figure 0. Spectra at different sites and for different 
field components show the same peaks, but with differing 
heights. Once the normal mode frequencies ujn are identi- 
fied, maps of the normal modes 5n{x) may be constructed 
by performing discrete Fourier sums on each component 
of the field as it evolves. For degenerate modes, each 
set of perturbed initial conditions gives a different lin- 
ear combination ti5i{x) of modes with the same fre- 
quency. Other linear combinations may also be produced 
by applying the symmetries of the static soliton 81,82, ■■■ 
to a given mode S. The degeneracy of a given frequency is 
found by computing the rank of the matrix of inner prod- 
ucts between the different linear combinations of degener- 
ate modes produced in these ways. Once the degeneracy 
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is determined, a complete orthonormal basis of modes 
with this frequency can be constructed. The character of 
any Oh group element can then be computed as a trace. 
These characters were within ±.001 of an integer value, 
and interpretation was unambiguous. 

Each peak in Figure |l| marks the frequency of a nor- 
mal mode. The lowest peak (at uj = .07) is the rotational 
zero mode, shifted to nonzero frequency by finite size ef- 
fects, effectively through interactions with image solitons 
one box length away. There are also several peaks cor- 
responding to relatively delocalised modes which we in- 
terpret as pion radiation. Two (at u — .545, .587) corre- 
spond to the lowest radiation modes, homogeneous away 
from the soliton, whose threefold degeneracy is split into 
_E+ + A2 by the presence of the soliton. The first radia- 
tion mode at nonzero wavenumber is at w = .908. The re- 
maining modes are the true vibrational excitations of the 
a particle. Somewhat fortuitously, the box size was small 
enough that the lowest inhomogeneous radiation mode 
has a frequency above the highest vibrational mode. 

Four widely separated B = 1 Skrymions have 24 zero 
modes, corresponding to 3 translations and 3 isorotations 
each. As they combine to form the B=4 soliton, 9 of them 
(3 translations, 3 rotations and 3 isorotations) remain as 
zero modes of the new system. One might expect the 
remaining 15 modes would survive as low energy vibra- 
tional modes. This is one fewer than what we find: we 
have an additional breather mode. 

The vibrational modes distort the B — 4 soliton as il- 
lustrated in Figure |^, and explained in the Table. The 
modes naturally divide into two sets. The lower nine vi- 
brational modes consist of deformations which, roughly, 
involve incompressible flow of the baryon charge. In 
contrast, the higher seven vibrational modes all have a 
'breathing' character, in which local baryon charge ex- 
pands or contracts to occupy a greater or lesser volume. 
The breather itself is simply a rescaling of the size of 
the soliton, with consequent change in density. The next 
mode up involves breathing motion of a dipole character, 
and the one above that of a quadrupolc nature. 

Remarkably, the vibrational modes below the breather 
fall into representations corresponding to those for small 
zero-mode deformations of the BPS 4-monopole solution 
[^. The same phenomenon occurs in the deuteron [p^ , 
and it will be straightforward to check the B = 3, 5, 6, 7 
solitons using the methods described here. Qualitative 
similarities between Skyrme multisolitons and BPS mul- 
timonopoles and their scattering dynamics have been 
noted before Our finding suggests a connection 

between the lowest energy Skyrmion vibrational modes 
and the multimonopole moduli spaces. If it holds for 
higher nuclei, there will be 45 — 7 such modes. It would 
be very interesting to interpret this number in terms of 
individual nucleon degrees of freedom (presumably trans- 
lations and spin/isospin). 

We have investigated the effect of doubling the param- 



eter on the spectrum of modes. All modes move up in 
frequency. The nine lowest modes move up by 15-25 per 
cent, the breathing modes by 30-45 per cent, and the ho- 
mogeneous pion modes roughly double in frequency. This 
has important phenomenological consequences. The low- 
est lying vibrational level for the real a particle is at 20.1 
MeV, whereas for — .526 our lowest lying frequency 
is at 0.7a;7r, or 94 MeV. For — 1.052 this is improved 
to 63 MeV. Whilst finite volume and grid spacing effects 
are present at a level of a few per cent in these results, we 
can safely conclude that a fit between the Skyrme model 
and nuclear gamma ray spectra will require parameters 
different than those usually used. It is interesting that 
a study of a single Skyrmion's breather mode, attempt- 
ing to identify it with the 'Roper resonance', also found 
that large values of ~ 1-5 were required ||l^. It re- 
mains to be seen whether such large values of can be 
accomodated in the theory. 
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Static B = 4 soliton, <^st (t>st +, - <5.367 4>st +, - <5|67 




FIG. 2. Contours of constant baryon density for the B=4 soliton, combined with its normal modes, as indexed by their 
frequencies in Figure |l|. These modes were studied in a box of size 8x8x8, with a grid spacing A^; = i in Skyrme units. For 
comparison, the soliton itself is a cube roughly 2x2x2 in these units. For the case uj = .367 two orthogonal modes of the 
degenerate multiplet are shown, in other cases a single mode only is shown. 



Frequency 


Degeneracy 


Symmetry 


Description 


.07 


3 




Rotations of the soliton. This is a zero mode broken by the finite box size. 


.367 


2 


E+ 


Lowest vibrational modes. One mode, 5^, alternately pulls the B — 4 cube into 
two 5 = 2 donuts in two perpendicular directions. The orthogonal mode, , 
pulls it into four B = 1 edges one way and two B = 2 donuts the other. 


.405 


1 




The corners of the cube make two interlacing tetrahedra. This mode pulls one 
tetrahedron out into four B = 1 corners, pushing the other one in. 


.419 


3 




Deform two opposite faces of the cube into rhombuses. 


.513 


3 


F2 


Deform the cube by pulling two opposite edges on one face, and the two perpen- 
dicular edges on the opposite face. This takes the cube to four B — 1 edges. 


.545 


2 


E+ 


Two of the pion fc = modes. 


.587 


1 




The remaining pion fc = mode, with tetrahedral symmetry. 


.605 


1 


At 


The breathing mode, with the full cubic symmetry of the soliton. 


.655 


3 


F2 


One face of the cube inflates, while the opposite face deflates. 


.738 


3 




One pair of diagonally opposite edges inflates, the parallel pair deflates. 


.908 


3 




The lowest nonzero {k — 1, 0, 0) pion radiative mode. 



TABLE L Description of the Skyrme fleld B = 4 normal modes marked in Figure 1. The notation of Hamermesh ||t| is used 
for the representations of the cubic group Oh', superscripts denote parity. 
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